Amplification of evanescent waves in a lossy left-handed material slab 
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We carry out finite-difference time-domain (FDTD) simulations, with a specially-designed bound- 
ary condition, on pure evanescent waves interacting with a lossy left-handed material (LHM) slab. 
Our results provide the first full-wave numerical evidence for the amplification of evanescent waves 
inside a LHM slab of finite absorption. The amplification is due to the interactions between the 
evanescent waves and the coupled surface polaritons at the two surfaces of the LHM slab and the 
physical process can be described by a simple model. 
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Recently, Pendry 1] suggested that a class of super- 
lenses could be realized by a planar slab of the so-called 
left-handed materials (LHMs) possessing both negative 
permittivity and permeability. With such a superlens, 
not only the phases of propagating waves but also the 
amplitudes of evanescent waves from a near-field object 
could be restored at its image. Therefore, the spatial 
resolution of the superlens could overcome the diffraction 
limit of conventional imaging systems and reach subwave- 
length scale. Especially at the single frequency where 
e(u>) = h(uj) = —1, the superlens could even be perfect 
- all the structural details of the object could be recon- 
structed at the image with unlimited resolution. It was 
soon realized that the superlens could be of great impor- 
tance in different technological areas. 

While great research interests were initiated by the rev- 
olutionary concept, hot debates were also raised 
5]. For example, based on the fact that LHMs are neces- 
sarily dispersive and dissipative 6, 7|, Garcia and Nieto- 
Vesperinas ||| argued that the unavoidable absorption in 
LHM, no matter how tiny it is, would suppress the am- 
plification of evanescent waves and thus the perfect lens 
is only a theoretical artifact and does not exist in real 
world. Despite the grea t amount of analytic and numer- 
ical attempts [IIIMEIEHQIIIO]!, reported 
results are diverse and disputes still remain. 

It is the purpose of this paper to examine the feasi- 
bility of the superlensing effect of a realistic LHM slab 
using a full-wave finite-difference time-domain (FDTD) 
method. Since it is the evanescent waves responsible for 
the superlensing effect, our investigation concentrates on 
the interaction between evanescent waves and the LHM 
slab. Our numerical results provide, for the first time, 
direct evidence that evanescent waves could be amplified 
inside a LHM slab of finite absorption. We show how the 
amplification happens and what is the role of the sur- 
face polaritons, excited by evanescent waves at the two 



surfaces of the LHM slab. We also study the effects of 
absorption and the thickness of the LHM slab on the 
amplification. 

The major innovation introduced in this paper is the 
incorporation of a specially-designed boundary condition 
in the FDTD scheme to simulate pure evanescent waves 
interacting with a LHM slab. We are therefore able to ex- 
plicitly study the characteristics of the evanescent waves 
and their dependence on different parameters. This is 
significantly different from previous works H Q, Hif . 
where both propagating and evanescent waves were in- 
volved in the simulation and the behavior of evanescent 
waves could only be implied by the image obtained, which 
caused large discrepancies in the conclusions since many 
other factors might have also contributed to the image 
quality in the FDTD simulations. 

In this paper we consider a two-dimensional transverse 
electric (TE) case. A LHM slab in vaccum extends from 
x = to L with surfaces normal to the rc-axis and is 
infinite in the yz-plane. The LHM is modelled by causal 
permittivity (e) and permeability (fi) with the identical 
plasmonic form, 



e(uj) = f/,(u>) = 1 — 



(1) 



where u> p is the plasma frequency and v c is the collision 
frequency. The z-polarized electric field source is excited 
in the x = — d (d = L/2 in our simulations) plane. The 
plane wave solution in vaccum is of the form, 



E z (x,y,t) = E z0 e 



— i(k x x-\-k v y—Ljt) 



(2) 



where k x and k y are wave numbers along the x- and y- 
directions, respectively, and they satisfy the dispersion 
relation, 



k 2 - — 



(3) 
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In the case of k 2 > k% (fcg = u 2 /c 2 ), k x is imaginary and 
corresponds to an evanescent wave along the x-direction. 

Experimentally, pure evanescent waves can be gener- 
ated by a prism using total internal reflection or a grating 
with a period finer than the wavelength of the radiation 
|l7|. Numerically, this can be done more simply in the 
FDTD simulation with the help of a specially-designed 
boundary condition along the transverse direction. Ow- 
ing to the continuity of the field across the interfaces, 
the transverse wave number k y is a constant in the whole 
space, both inside and outside the LHM slab. Therefore, 
it is possible for us to explicitly assign a predefined k y for 
the simulated wave. One way to achieve this objective 
is to impose a boundary condition in y-direction. In this 
case, we require that the fields in the boundary satisfy 
E z (x,y ± Ay) = E z {x, y)e TlkyAy . This approach is an 
extension of the periodic boundary condition fPBC")[Isj 
where, normally, only the wave vectors within the first 
Brillouin zone are employed. 

In the FDTD simulation, we apply the standard Yee's 
algorithm in which the leapfrog staggered electric field 
and magnetic field are updated in time using the two cou- 
pled Maxwell's curl equations [l^. The computational 
space is taken to be 4000 x 1 cells with Ax = Ay = 0.3 
mm. At both ends in the x-direction, 10-cell uniax- 
ial anisotropic perfectly-matched layer (UPML) ^| is 
added to truncate the computational space. The time 
step used is At = Axj '(2c) = 0.5 ps. A sinusoidal source 
with a smooth turn-on of 30 periods Q is employed to 
excite a monochromatic wave at the angular frequency 
u> w 11 GHz. The corresponding wavelength in vacuum 
Ao (=27r/fco) is about 566Ax. By choosing lo p and v c in 
(1) appropriately, we have e(u>) — /J>(u>) = —I — ij (7 de- 
notes the loss term) in all our simulations. The dispersive 
e and /i is handled in time domain using the piecewise- 
linear recursive convolution (PLRC) method (lijj . 

In the simulation, we monitored the time evolution of 
E z at three different locations, the two surfaces of the 
LHM slab and the image plane. It was found that a 
transient process occurs at the begining of the simula- 
tion and after a certain period, the fields reach their 
steady state and vibrate sinusoidally with unchanged am- 
plitudes. The time required to reach the final steady 
state and the dynamic characteristics of the time evolu- 
tion vary dramatically with different parameters (7, k y 
and L). It is beyond the scope of this paper to describe 
the details of these simulations. Instead, we will only 
discuss the steady-state results obtained. 

Figure 1(a) shows the spatial distribution of the 
steady-state amplitude of E z as a result of an evanes- 
cent wave (k y = 2fco) interacting with a LHM slab of 
thickness L = 80 Ax ~ 0.14Ao for different values of 7, 
ranging from 0.001 to 0.2. Figure 1(b) shows the corre- 
sponding E z amplitude at the two surfaces of the slab 
as a function of 7. The most important feature revealed 
by the figures is the amplification of the evanescent wave 
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FIG. 1: An evanescent wave (k y = 2ko) interacting with a 
LHM slab (L = 80 Ax ~ 0.14Ao) for different values of absorp- 
tion (7), ranging from 0.001 to 0.2. (a) The x-distribution of 
the steady-state E z amplitude; (b) the steady-state E z am- 
plitude at the two surfaces of the LHM slab as a function of 
7. The vertical thin dot-lines in (a) denote the two surfaces 
of the LHM slab. 



within the LHM slab, i.e., the field at the second surface 
is stronger than that at the first surface, for finite ab- 
sorption, even though absorption tends to suppress the 
amplification. For small values of 7 from 0.001 to 0.01, 
the field distribution is almost identical. The evanescent 
wave from the source undergos very little reflection at 
the first surface of the LHM slab and the field between 
the source and the first surface is only slightly different 
from the exponential decay expected. Meanwhile, the 
field inside the slab shows a monotonic increase which is 
close to the exponential increase predicted by Pendry for 
the ideal case (e = fi = —1) 1]. The field is amplified 
by about 4 times from 0.55 at the first surface to 2.1 at 
the second surface. When 7 increases above 0.01, the 
field at the first surface increases noticeably while that 
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FIG. 2: An evanescent wave (k y = 2ko) interacting with a 
LHM slab (L = 160Az » 0.28A ) for different values of ab- 
sorption (7), ranging from 0.001 to 0.2. (a) The ^-distribution 
of the steady-state E z amplitude; (b) the steady-state E z am- 
plitude at the two surfaces of the LHM slab as a function of 
7. The vertical thin dot-lines in (a) denote the two surfaces 
of the LHM slab. 



at the second surface drops, so that the amplification is 
weakened. When 7 reaches 0.2, the amplification still 
exists but now the field at the second surface (=1.56) is 
only about 1.5 times as strong as that of the first surface 
(=0.97). From the trends of the curves shown in Figure 
1(b), we can expect that if 7 is further increased above a 
certain value 7 C (> 0.2), the field at the first surface may 
finally increase above the decreasing field at the second 
surface, so that the amplification would be replaced by 
decay. Associated with the increase in the field at the 
first surface, larger reflections of the evanescent wave oc- 
cur at the first surface. Furthermore, a peak appears at 
the first surface for 7 > 0.1 and the peak becomes more 
prominent with increasing 7. 

Figure 2 shows the simulation results for a thicker slab 



with L — 160Ax ~ 0.28Ao while keeping k y = 2ko and 
7 from 0.001 to 0.2. While the results show some com- 
mon features with those obtained for the thinner slab 
presented in Figure 1, they also reveals a number of dis- 
tinct features. For 7 = 0.001, the LHM slab amplifies 
the field from 1.03 at the first surface to 4.06 at the sec- 
ond surface. The amplification is gradually weakened 
with increasing field at the first surface and decreasing 
field at the second surface while 7 is raised. At around 
7 = 0.07 (7c), a crossover occurs where the fields at the 
two surfaces become equal. For 7 > 0.07, the field at the 
first surface is stronger than that at the second surface, 
and the amplification turns into decay. Compared with 
7 C > 0.2 obtained for the thinner slab, the crossover ab- 
sorption value 7c = 0.07 is lowered dramatically for the 
thicker slab. The smaller 7 C obtained for the thicker slab 
suggests that large L is not favored for the amplification 
of the evanescent wave inside the LHM slab since it puts 
more stringent constraints. It is also clearly seen from 
Figure 2(a) that a peak is present, and large reflection 
occurs at the first surface, even for the lowest absorption 
studied (7 = 0.001) where the LHM slab is supposed to 
be well matched to vacuum. This is different from the 
results shown in Figure 1(a). 

So far we have shown the simulation results of the 
evanescent wave interacting with LHM slabs of differ- 
ent absorption and thickness. The numerical examples 
provide direct evidence that an evanescent wave could 
be amplified in a LHM slab of finite absorption. As a 
result of the amplification of evanescent waves inside the 
lossy LHM slab, superlensing effect with subwavelcngth 
image resolution could be achieved practically |2fJ. It 
should be emphasized that our results are based on full- 
wave numerical experiments where only the casual e(w) 
and h(<jj) of the LHM are predefined and the results are 
obtained without any further assumptions and approxi- 
mations. 

Next, we will discuss the physical origin of the ampli- 
fication of the evanescent waves inside the LHM slabs. 
It was shown by Ruppin that a vacuum/LHM interface 
supports surface polaritons (SPs), which are character- 
ized by exponential decay in field on both sides of the 
interface |2lj, |22j ■ The peaks at the two surfaces of the 
LHM slabs and their shapes, shown in Figure 1(a) and 
2(a), clearly illustrate the presence of SPs and emphasize 
their important roles in the systems. We believe that the 
coupled SPs, excited by the evanescent waves at the two 
surfaces of the LHM slab, are responsible for the ampli- 
fication observed 0, 0, To prove this, we show 
that the features revealed by our numerical results are 
consistent with the predictions of a simple model of the 
interaction between the evanescent wave and the coupled 
SPs. 

The model we use to describe the interaction between 
the evanescent wave and the two coupled SPs is simply 
the same as the forced vibration and resonance of damped 
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identical oscillators with linear coupling;, which can be 
written in the coupled equations |l3| 



701 

'702 



Wq01 
'^002 



n 2 c 



Fe iuJt 




(4) 



where (pi and <p2 denote the oscillations of the first SP 
and the second SP, respectively; ujq is the natural reso- 
nant frequency of the single SP; tt 2 describes the coupling 
strength between the two SPs; and F describes the exter- 
nal excitation introduced by the evanescent wave. When 
the evanescent wave is of frequency loq, which is the case 
in our simulation, the steady-state solutions of <p\ and (f>2 
have the simple forms 



4>i (^o) 



-ijoj Fe luJot 
^ + 7 2 ^o 2 ' 



4>2 (Wo) 



Q 2 Fe 



(5) 



The characteristics of the coupled SPs and their depen- 
dences on different parameters can be readily analyzed 
from these solutions. In the lossless limit 7 — > 0, the first 
SP disappears with |0i| = and only the second SP is 
present in the system with finite coupling. This is the 
ideal case and an exponential amplification occurs inside 
the LHM slab. For finite absorption both SPs are ex- 
cited. Simple calculation shows that, \4>i \ increases with 
increasing 7 until it passes through a peak, located at 
7 = 7 P (jp = f^/wo), and decreases beyond. On the 
other hand, \4>2\ decreases monotonically with increasing 
7. The ratio \<j>i/<j>2\ — 7^o/^e gives a measure of the 
relative strength of the two SPs. For small values of 7, 
|0i/02| < 1, the second SP is stronger than the first so 
that amplification takes place. However, the effect is re- 
versed for 7 > 7c (7c = ^c/ w o) where the first SP is 
dominant. It is clear that 7 C = 7p, which means that the 
crossover of \<fii\ and \<p2\ is accompanied by the peak in 
\cf)i\. All these features are clearly observed in Figures 1 
and 2. Furthermore, it is expected that increasing L re- 
duces Q 2 , so that the value of 7 C is smaller for the thicker 
LHM slab, which is also consistent with the numerical re- 
sults. 

To explain why the first SP is stronger for the thicker 
slab when 7 is small and unchanged, we need to know the 
expressions of fl 2 and F explicitly on the geometry. In 

the weak coupling limit k x L 



00 (k x = J, 



k 2 — U! 2 /c 2 ), 



they can be approximately written in the forms 

urn 

F fti Ce~ KlL l 2 and £l 2 c « De~ KxL , where C and D are 
coefficients independent of k x L. Then, we have 



d|<M^o)| 
8{kL) 



1 



JLUqF 



2(04 + 7^ 



2 2\ 

^o7 ) 



(6) 



When 30* - u?^ 2 > 0, i.e., 7 < a/^c, the first SP is 
getting stronger with increasing L, which explains the 
stronger first SP presented in Figure 2. It should be men- 
tioned that the weak coupling limit may not be applicable 
for our simulations, where k x L ps 0.24 (for L = 8OA2;) 



and 0.48 (for L = I6OA2;). However, it aids in un- 
derstanding the result. In fact, from a purely physical 
point of view, the emergence of the first SP with certain 
strength is crucial for the effective excitation of the sec- 
ond SP in the thicker slab, where the coupling coefficient 
VL 2 is reduced. It is the first SP that mediates the energy 
transfer from the evanescnet waves to the second SP. 

In conclusion, we have proposed a method to simulate 
pure evanescent waves interacting with a LHM slab using 
FDTD. The simulation results provide the first full- wave 
numerical evidence that evanescent waves could be am- 
plified inside a LHM slab with finite absorption. It also 
reveals that the coupled SPs, excited by the evanescent 
waves at the two surfaces of the LHM slab, play impor- 
tant roles in the amplification effect. We have shown that 
absorption tends to suppress the amplification and the 
amplification could be replaced by decay when the ab- 
sorption is above the crossover value 7 C , which is greatly 
reduced with increasing thickness of the LHM slab. All 
features in our numerical results can be understood con- 
sistently using a simple model which describes the in- 
teraction between the evanescent waves and the coupled 
SPs. 
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